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ABSTRACT 

Departures  from  axisymmetric  balance  dynamics  are  quantified  during  a  case  of  secondary  eyewall  for¬ 
mation.  The  case  occurred  in  a  three-dimensional  mesoscale  convection-permitting  numerical  simulation  of 
a  tropical  cyclone,  integrated  from  an  initial  weak  mesoscale  vortex  in  an  idealized  quiescent  environment. 
The  simulation  exhibits  a  canonical  eyewall  replacement  cycle.  Departures  from  balance  dynamics  are 
quantified  by  comparing  the  azimuthally  averaged  secondary  circulation  and  corresponding  tangential  wind 
tendencies  of  the  mesoscale  integration  with  those  diagnosed  as  the  axisymmetric  balanced  response  of 
a  vortex  subject  to  diabatic  and  tangential  momentum  forcing.  Balance  dynamics  is  defined  here,  following  the 
tropical  cyclone  literature,  as  those  processes  that  maintain  a  vortex  in  axisymmetric  thermal  wind  balance. 

The  dynamical  and  thermodynamical  fields  needed  to  characterize  the  background  vortex  for  the  Sawyer- 
Eliassen  inversion  are  obtained  by  azimuthally  averaging  the  relevant  quantities  in  the  mesoscale  integration 
and  by  computing  their  corresponding  balanced  fields.  Substantial  differences  between  azimuthal  averages 
and  their  homologous  balance-derived  fields  are  found  in  the  boundary  layer.  These  differences  illustrate  the 
inappropriateness  of  the  balance  assumption  in  this  region  of  the  vortex  (where  the  secondary  eyewall  tan¬ 
gential  wind  maximum  emerges).  Although  the  balance  model  does  broadly  capture  the  sense  of  the  forced 
transverse  (overturning)  circulation,  the  balance  model  is  shown  to  significantly  underestimate  the  inflow  in 
the  boundary  layer.  This  difference  translates  to  unexpected  qualitative  differences  in  the  tangential  wind 
tendency.  The  main  finding  is  that  balance  dynamics  does  not  capture  the  tangential  wind  spinup  during  the 
simulated  secondary  eyewall  formation  event. 


1.  Introduction 

A  secondary  eyewall  of  a  tropical  cyclone  is  a  prom¬ 
inent  coherent  structure  that  is  concentric  to  the  primary 
eyewall  and  is  characterized  by  maxima  in  convective 
activity  and  tangential  winds.  Secondary  eyewalls  are 
common  in  the  most  intense  tropical  cyclones  (Hawkins 
and  Helveston  2004,  2008;  Kossin  and  Sitkowski  2009; 
Kuo  et  al.  2008).  Secondary  eyewalls  have  been  found  to 
be  associated  with  intensity  changes  (Yang  et  al.  2013), 
including  eyewall  replacement  cycles  (Houze  et  al.  2007, 
and  references  therein),  and  with  the  longer  duration  of 
higher  storm  intensity  (Kuo  et  al.  2008)  and  storm 
growth  (Maclay  et  al.  2008).  Despite  the  recognized 
importance  of  secondary  eyewalls,  there  is,  as  of  yet,  no 
unified  theory  to  explain  secondary  eyewall  formation 
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(SEF)  and  to  guide  the  development  of  forecast  tools, 
which  today  tend  to  rely  on  empirical  relations  (Kossin 
and  Sitkowski  2009). 

A  wide  variety  of  physical  processes  aiming  to  explain 
SEF  have  been  described.  These  include  two-dimensional 
vortex  interactions  (Kuo  et  al.  2004,  2008),  anisotropic 
upscale  energy  cascades  (Terwey  and  Montgomery  2008), 
and  changes  in  axisymmetric  efficiency  [as  defined  by 
Schubert  and  Hack  (1982)  and  Hack  and  Schubert  (1986)] 
related  to  the  radial  expansion  of  the  azimuthally  aver¬ 
aged  tangential  wind  field  and  sustained  latent  heating 
outside  of  the  primary  eyewall  (Rozoff  et  al.  2012).  Other 
models  propose  wave-mean  flow  interaction  associated 
primarily  with  vortex  Rossby  waves  and  their  modifi¬ 
cation  to  the  storm  inner-core  region  (Montgomery  and 
Kallenbach  1997;  Martinez  et  al.  2010a, b;  Menelaou 
et  al.  2012;  Qiu  et  al.  2010;  Abarca  and  Corbosiero  2011) 
or  the  generation  and  accumulation  of  convectively  in¬ 
duced  potential  vorticity  anomalies  and  the  absence  of 
vortex  Rossby  wave  activity  (Judt  and  Chen  2010).  Sun 
et  al.  (2013)  proposed  that  SEF  arises  from  a  feedback 
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between  inward-moving  rainbands,  the  balanced  re¬ 
sponse  to  their  source  of  heat,  and  unbalanced  dynam¬ 
ics.  Wang  et  al.  (2013)  proposed  that  a  “net”  radial  force 
in  the  boundary  layer  is  relevant  in  the  process  and  Qiu 
and  Tan  (2013)  emphasized  the  role  of  the  asymmetric 
boundary  layer.  Kepert  (2013)  presented  theoretical 
arguments,  based  mostly  on  linearized  Ekman  theory, 
to  propose  a  feedback  mechanism  for  secondary  eye- 
wall  formation  that  involves  a  local  enhancement  of  the 
radial  vorticity  gradient,  frictional  updraft,  and  con¬ 
vection.  The  hypothesized  feedback,  based  primarily 
on  linear  Ekman  balance  reasoning,  has  been  falsified 
by  Montgomery  et  al.  (2014). 

SEF  has  been  proposed  recently  as  an  element  of 
a  new  paradigm  of  tropical  cyclone  intensification 
(Huang  et  al.  2012;  Abarca  and  Montgomery  2013).  As 
part  of  the  new  paradigm,  Smith  et  al.  (2009)  showed 
that  from  the  azimuthally  averaged  perspective,  the 
spinup  of  the  system-scale  swirling  circulation  (i.e.,  the 
radial  convergence  of  absolute  angular  momentum) 
occurs  via  two  mechanisms.  The  first  mechanism  is 
physically  rooted  above  the  boundary  layer,  in  the  sense 
that  it  is  the  result  of  the  radial  gradient  of  buoyancy 
generation  rate,  associated  with  the  aggregate  latent 
heat  release  in  deep  convective  clouds.  This  mechanism 
can  be  captured  approximately  by  balance  dynamics  [as 
defined  by  Shapiro  and  Willoughby  (1982)]  and  has  been 
formulated  in  the  context  of  the  conservation  of  abso¬ 
lute  angular  momentum.  The  second  spinup  mechanism 
is  physically  rooted  within  the  boundary  layer  and  is  the 
result  of  inward  radial  advection  of  absolute  angular 
momentum  (associated  with  frictional  inflow)  surpass¬ 
ing  the  depletion  of  absolute  angular  momentum  by 
frictional  torque.  While  this  latter  mechanism  cannot 
operate  without  the  first  mechanism  (through  the  cou¬ 
pling  with  the  interior  flow  via  the  radial  pressure  gradi¬ 
ent  force),  it  is  essentially  a  mechanism  of  the  dynamics  of 
the  boundary  layer  and  is  manifested  as  a  type  of  vortex 
breakdown.  Accordingly,  the  boundary  layer  spinup 
mechanism  lies  outside  the  realm  of  balance  dynamics.  In 
this  study,  as  in  Abarca  and  Montgomery  (2013),  we 
follow  Smith  et  al.  (2009)  and  Smith  and  Montgomery 
(2010)  in  the  adoption  of  a  dynamical  definition  of  the 
boundary  layer.  The  term  boundary  layer  will  be  used 
here  to  describe  the  shallow  layer  of  strong  inflow  near 
the  sea  surface  that  arises  largely  because  of  the  frictional 
disruption  of  balance  dynamics  near  the  surface  (e.g., 
Smith  et  al.  2009,  their  Fig.  6). 

Scale  analysis,  high-resolution  mesoscale  integrations, 
and  in  situ  observations  of  the  hurricane  (Willoughby 
1979;  Smith  et  al.  2009;  Moller  and  Shapiro  2002;  Zhang 
et  al.  2001;  Willoughby  1990;  Bell  and  Montgomery 
2008)  strongly  suggest  that  on  the  vortex  scale,  a  large 


portion  of  the  hurricane  flow  is  in  axisymmetric  gradient 
and  hydrostatic  balance  and  consequently  this  part  of 
the  vortex  should  satisfy  the  thermal  wind  equation.  To 
stay  in  such  a  balance  state,  a  secondary  circulation 
will  develop  to  oppose  the  forcing  of  diabatic  heating 
and/or  friction.  After  the  seminal  work  of  Eliassen 
(1951),  a  number  of  authors  (e.g.,  Willoughby  1979; 
Shapiro  and  Willoughby  1982)  have  used  the  so-called 
Sawyer-Eliassen  equation  to  diagnose  the  secondary 
circulation  of  a  hurricane-like  vortex  forced  by  diabatic 
heating  or  momentum  sources.  The  Sawyer-Eliassen 
equation  has  been  used  also  as  a  mean-flow  model  in 
the  limit  of  weak  departures  from  circular  flow,  to  study 
axisymmetric  tropical  cyclone  intensification  (e.g.,  Bui 
et  al.  2009),  tropical  cyclone  asymmetries  (e.g.,  Shapiro 
and  Montgomery  1993;  Montgomery  and  Shapiro  1995; 
Moller  and  Montgomery  2000;  Moller  and  Shapiro 
2002),  the  interaction  of  the  tropical  cyclone  with  its 
environment  (e.g.,  Challa  and  Pfeffer  1980;  Pfeffer  and 
Challa  1981;  Persing  et  al.  2002;  Moller  and  Shapiro 
2002),  and  the  extent  to  which  vortex  evolution  during 
tropical  cyclogenesis  proceeds  via  axisymmetric  balance 
dynamics  (Hendricks  et  al.  2004;  Montgomery  et  al. 
2006;  Fang  and  Zhang  2011;  Wang  2012). 

As  stated  in  the  foregoing  discussion,  the  azimuthally 
averaged  flow  associated  with  the  first  mechanism  for 
system-scale  spinup  can  be  described  approximately  by 
balance  dynamics.  The  more  general  conception  of  bal¬ 
ance  flow  is  a  flow  in  which  the  velocity  field  is  func¬ 
tionally  related  to  the  mass  field  (presumed  to  be  related 
quasi-statically  to  the  pressure  field;  McIntyre  2012). 
However,  in  this  study,  following  the  custom  in  tropical 
cyclone  research,  axisymmetric  balance  flow  will  be  de¬ 
fined  as  that  part  of  the  hurricane  primary  circulation 
that  satisfies  the  gradient  and  hydrostatic  balance.  We 
refer  to  axisymmetric  balance  dynamics  as  the  axisym¬ 
metric  processes  that  constrain  the  flow  to  remain  in  axi¬ 
symmetric  thermal  wind  balance.  Here,  the  axisymmetric 
portion  of  the  hurricane  flow  that  does  not  satisfy  the 
thermal  wind  equation  is  called  unbalanced  flow  and  the 
dynamics  associated  with  unbalanced  flow  are  referred 
to  as  unbalanced  dynamics.  Unlike  many  other  studies, 
the  concept  of  unbalance  does  not  refer  primarily  to  the 
effect  of  gravity  waves.  Instead,  it  includes  the  strong 
horizontal  advective  dynamics  that  arise  because  the 
momentum  equations  are  coupled,  on  account  of  fric¬ 
tion,  in  the  region  spanning  the  transition  layer  between 
the  surface  layer  and  the  nearly  inviscid  vortex  circula¬ 
tion  in  the  bulk  flow  aloft. 

It  has  been  proposed  that  unbalanced  dynamics,  as 
defined  herein,  are  an  important  element  in  SEF.  Using 
a  model-observation  consistent  dataset,  Wu  et  al.  (2012) 
and  Huang  et  al.  (2012)  proposed  that  SEF  is  the 
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culmination  of  a  sequence  of  structural  changes  of  the 
mature  tropical  cyclone’s  inner  core.  These  changes  start 
with  a  radial  expansion  of  the  tangential  wind  field,  fol¬ 
lowed  by  the  response  of  the  boundary  layer,  character¬ 
ized  by  an  increase  in  radial  inflow  and  subsequently  by 
the  development  of  supergradient  winds  and  an  increase 
in  horizontal  convergence.  Abarca  and  Montgomery 
(2013)  presented  a  verification  of  this  hypothesis  using 
independent  evidence  from  the  idealized  full-physics, 
high-resolution  integration  originally  studied  by  Terwey 
and  Montgomery  (2008)  and  by  Terwey  et  al.  (2013)  and 
with  a  slab  boundary  layer  model.  With  this  evidence, 
Abarca  and  Montgomery  (2013)  proposed  that  un¬ 
balanced  dynamics  are  an  important  element  for  SEF. 
Alongside  the  model  evidence,  there  is  also  in  situ  ob¬ 
servational  evidence  that  mature  secondary  eyewalls  ex¬ 
hibit  supergradient  winds  near  the  top  of  the  boundary 
layer  (Didlake  and  Houze  2011;  Bell  et  al.  2012). 

In  this  study  we  will  present  a  quantification  of  the 
degree  to  which  SEF  departs  from  balanced  dynamics.  It 
is  based  on  the  canonical  SEF  integration  studied  by 
Terwey  and  Montgomery  (2008),  Terwey  et  al.  (2013), 
and  Abarca  and  Montgomery  (2013). 

The  manuscript  is  organized  as  follows:  section  2 
describes  the  scientific  methodology,  the  models,  and 
their  integrations;  section  3  presents  the  main  results, 
and  section  4  provides  a  summary  of  the  findings  and 
concluding  remarks.  Two  appendixes  are  included  to 
address  important  methodological  questions  regarding 
the  balance  integrations  and  their  robustness  to  the 
characterization  of  the  background  vortex. 

2.  Scientific  methodology,  models,  and  integrations 

The  methodology  of  this  work  consists  essentially  of 
comparing  the  results  of  two  different  flow  model  frame¬ 
works.  One  modeling  framework  is  the  Regional  Atmo¬ 
spheric  Modeling  System  (RAMS;  Pielke  et  al.  1992; 
Cotton  et  al.  2002).  The  other  is  the  Sawyer-Eliassen 
equation  and  corresponding  balanced  tendency  pre¬ 
diction  as  presented  in  Bui  et  al.  (2009).  We  use  azimuthal 
averages  of  the  high-resolution  model  output  of  a  ca¬ 
nonical  SEF  event  to  integrate  the  Sawyer-Eliassen 
equation.  A  comparison  of  the  secondary  circulation 
and  the  derived  tangential  wind  tendency  of  the  two 
modeling  frameworks  is  performed  to  quantify  the 
extent  to  which  the  high-resolution,  full-physics  RAMS 
integration  undergoing  SEF  is  captured  by  axisym- 
metric  balance  dynamics.  This  procedure  has  been  used 
previously  to  study  the  formation  stage  of  incipient  hur¬ 
ricane  vortices  (e.g.,  Hendricks  et  al.  2004;  Montgomery 
et  al.  2006;  Wang  2012)  and  in  mature  hurricanes  (e.g., 
Bui  et  al.  2009;  Fudeyasu  and  Wang  2011).  The  remainder 


of  this  section  summarizes  the  essential  features  of 
the  two  models  and  specific  details  of  the  methodology 
used  to  quantify  departures  from  axisymmetric  balance 
during  SEF. 

a.  The  RAMS  and  the  RAMS  integration 

RAMS  is  a  three-dimensional  nonhydrostatic  nu¬ 
merical  modeling  system.  For  the  purpose  of  the  present 
study,  we  summarize  the  numerical  integration,  which 
exhibits  a  canonical  SEF  and  eye  wall  replacement  cycle, 
presented  by  Terwey  and  Montgomery  (2008)  and  re¬ 
visited  by  Terwey  et  al.  (2013),  Abarca  and  Montgomery 
(2013),  and  Montgomery  et  al.  (2014).  This  simulation 
uses  a  surface  flux  parameterization  based  on  the  Louis 
(1979)  scheme;  the  radiation  scheme  introduced  by 
Harrington  (1997);  a  subgrid-scale  turbulence  scheme 
based  on  Smagorinsky  (1963),  with  the  modifica¬ 
tions  from  Lilly  (1962)  and  Hill  (1974);  and  the  seven- 
species  microphysical  scheme  by  Walko  et  al.  (1995),  with 
the  specification  described  in  Montgomery  et  al.  (2006). 

The  RAMS  simulation  used  in  this  study  is  an  idealized 
integration  on  an /plane  (15°N)  with  a  constant  sea  sur¬ 
face  temperature  of  28°C.  The  model  is  configured  with 
30  vertical  levels  (vertical  grid  spacing  varies  from  300  m 
near  the  surface  to  1800  m  near  the  top  of  the  domain,  at 
about  26-km  altitude).  The  model  configuration  has 
a  parent  domain,  and  two  (two-way  interactive)  nested 
domains.  Each  nested  domain  is  located  at  the  center  of 
the  parent  grid.  Horizontal  grid  spacing  is  24, 6,  and  2  km; 
and  there  are  168,  170,  and  251  grid  points  respectively. 
The  outermost  grid  is  cyclic  in  the  horizontal  direction. 
The  initial  condition  is  characterized  by  a  quiescent  en¬ 
vironment  [Jordan’s  (1958)  sounding].  It  contains  a  weak 
(10ms_1)  mesoscale  cyclonic  vortex,  located  at  the  cen¬ 
ter  of  the  three  domains  (with  maximum  speed  at  4-km 
height  and  75-km  radius),  that  is  initially  in  thermal  wind 
balance.  A  positive  moisture  anomaly  (water  vapor 
mixing  ratio  increased  up  to  1.3  g  kg-1  near  the  center  of 
the  initial  vortex)  is  specified  also.  Further  details  of  the 
RAMS  experimental  setup  can  be  found  in  Terwey  and 
Montgomery  (2008),  Montgomery  et  al.  (2006),  Terwey 
et  al.  (2013),  and  Abarca  and  Montgomery  (2013).  In  this 
study  we  present  azimuthal  and  1-h  time  averages  of  the 
simulated  variables.  Such  averages  are  obtained  after 
interpolating  the  RAMS  data  from  its  original  Arakawa 
C  grid  in  Cartesian  coordinates  into  a  cylindrical  grid  (the 
center  is  determined  as  the  centroid  of  the  potential 
vorticity  field  in  the  lowest  7.3  km  of  the  domain). 

b.  The  Sawyer-Eliassen  model 

The  Sawyer-Eliassen  equation  is  a  diagnostic  equation 
for  the  meridional  streamfunction  in  the  radius-height 
plane.  This  equation  determines  the  meridional  circulation 
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required  to  maintain  a  vortex,  under  heat  and  momen¬ 
tum  forcing,  in  hydrostatic  and  gradient  wind  balance. 
There  are  several  derivations  of  the  Sawyer-Eliassen 
equation  in  the  literature  (e.g.,  Shapiro  and  Willoughby 
1982;  Schubert  et  al.  2007).  In  this  study  we  integrate  the 
version  derived  by  Bui  et  al.  (2009): 


d 

dr 


dx  1  dip 
g  dz  pr  dr 


1  dip 
pr  dz_ 
1  dip 
pr  dz 


=  Syr(x2Q)  +  yz (Cx2Q)  -  / (xiFx)-  (1) 


Here  g  represents  the  acceleration  due  to  gravity 
(9.81  ms-2), /is  the  Coriolis  parameter  (evaluated  here 
at  15°N,  as  in  the  RAMS  simulation),  p  is  the  air  density, 
r  is  the  radius,  and  z  is  the  physical  height  above  sea 
level.  Also,  x  =  1/0,  where  6  is  the  potential  tempera¬ 
ture,  f  =  (1  /r)(drv/dr)  is  the  vertical  relative  vorticity  for 
the  mean  vortex  (i v  is  the  tangential  wind  velocity),  C  = 
(v2/r)  +  fu  is  the  compound  Coriolis  and  centrifugal  force 
per  unit  mass,  and  £  =  (2 v/r)  +  /is  the  modified  Coriolis 
parameter.  The  transverse  streamfunction  ip  is  defined 
such  that  the  radial  and  vertical  wind  velocities  u  and  w, 
respectively,  are  u  =  -(1  lr)(dip/dz)  and  w  =  (l/r)(dip/dr) 
[see  Bui  et  al.  (2009)  for  details].  Finally,  Q  and  Fx  are  the 
forcing  functions,  representing  the  diabatic  heating  and 
tangential  momentum  sources,  respectively.  These  last 
two  functions  are  defined  explicitly  below,  in  section  2c(2), 
along  with  details  of  the  integration  of  Eq.  (1). 


c.  The  Sawyer-Eliassen  integrations 

To  characterize  the  balanced  vortex  and  the  forcing 
functions  used  to  integrate  the  balance  model,  we  use 
RAMS  azimuthal  and  1-h  averaged  fields  centered  at  hour 
174  (when  the  simulation  was  undergoing  SEF;  see  section 
3  for  details  of  the  secondary  eyewall  in  the  mesoscale 
model).  In  this  section  we  describe  how  we  integrate  the 
Sawyer-Eliassen  model,  using  the  RAMS  integration.  We 
describe  also  aspects  regarding  the  ellipticity  condition  in 
the  Sawyer-Eliassen  integrations. 


thermal  wind  equation,  as  well  as  the  equation  of  state  to 
calculate,  without  approximation,  the  thermodynamical 
fields  necessary  to  integrate  Eq.  (1). 

The  second  method  of  constructing  the  balanced  vortex 
is  to  use  the  azimuthally  averaged  pressure  and  density 
fields  from  the  RAMS  model  output  to  infer  the  corre¬ 
sponding  balanced  tangential  wind  field.  In  this  method 
we  follow  the  definition  of  gradient  wind  flow  in  Holton 
(2004).  As  pointed  out  in  Bui  et  al.  (2009),  for  a  given 
mass  field,  it  is  not  always  possible  to  calculate  the  cor¬ 
responding  balanced  tangential  wind  field.  This  com¬ 
putation  requires  the  solution  of  a  quadratic  formula  in 
which  the  radicand  can  become  negative.  In  the  data 
used  in  this  work,  the  radicand  becomes  negative  only  in 
some  sparse  regions  of  the  outflow  layer  and  the  prob¬ 
lem  is  avoided  locally  by  setting  dp! dr  =  0  when  the 
quantity  becomes  positive. 

At  this  point,  it  is  useful  to  recall  that  azimuthally 
averaged  thermodynamical  and  dynamical  fields  from 
mature  hurricane  mesoscale  integrations  are  not  gen¬ 
erally  in  axisymmetric  balance  with  each  other.  How¬ 
ever,  it  has  been  a  common  practice  to  integrate  the 
Sawyer-Eliassen  model  with  such  azimuthal  averages 
(e.g.,  Hendricks  et  al.  2004;  Montgomery  et  al.  2006; 
Fang  and  Zhang  2011;  Wang  2012).  In  appendix  B  we 
present  an  integration  of  the  Sawyer-Eliassen  model 
characterizing  the  vortex  with  azimuthally  averaged 
fields  from  the  RAMS  integration.  It  will  be  shown  here 
that  the  conclusions  obtained  using  the  strict  balance 
definitions  of  Smith  (2006)  and  Holton  (2004))  or  azi¬ 
muthal  averages  are  qualitatively  consistent  with  each 
other  so  that  the  general  conclusions  of  this  manuscript 
are  robust  to  any  of  these  methodologies.1 

2)  Forcing  functions 

The  Sawyer-Eliassen  model  is  integrated  with  dia¬ 
batic  heating  and  momentum  forcing  diagnosed  from 
azimuthal  and  time  averages  of  the  RAMS  integrations. 
Following  Bui  et  al.  (2009),  the  azimuthally  averaged 
and  time-averaged  diabatic  forcing  Q  is  computed  as 
the  difference  between  the  total  diabatic  heating  rate 
and  eddy  heat  fluxes,  ignoring  subgrid-scale  heat 
diffusion,  as 


1)  Balanced  vortex 

The  vortex  used  to  integrate  the  balance  model  is 
constructed  here  in  two  different  ways  (and  other  ways  in 
appendix  B).  First,  we  follow  Smith  (2006)  and  use  the 
RAMS  azimuthally  averaged  and  time-averaged  tangen¬ 
tial  wind  field  to  compute  its  corresponding  balanced 
density  field.  In  this  procedure,  to  initialize  the  density 
integration,  we  use  the  vertical  density  profile  from  the 
RAMS  model  at  large  radius.  This  methodology  uses  the 


Q-i-**- 

v  dr 


v'  30'  ,30' 

- —  —  w — . 

r  dX  dz 


(2) 


1  Although  it  should  be  clear  from  the  findings  summarized 
herein,  we  point  out  that  the  robustness  of  our  findings  should  not 
be  construed  to  mean  that  the  azimuthally  averaged  fields  are  in 
thermal  wind  balance.  This  is  indeed  not  the  case,  as  the  differences 

displayed  in  Figs.  2  and  3  attest.  See  upcoming  discussion  for  fur¬ 
ther  details. 
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Here  the  overbar  denotes  azimuthal  average  on 
a  constant  height  surface  and  the  prime  denotes  de¬ 
viations  therefrom;  6  represents  the  azimuthally  aver¬ 
aged  total  diabatic  heating  rate.  Other  symbols  are  as 
defined  previously. 

The  azimuthally  averaged  and  time-averaged  tangen¬ 
tial  momentum  source  F  is  computed  as  the  difference 
between  the  local  time  change  of  tangential  momentum 
and  the  sum  of  the  mean  absolute  vorticity  flux  and  the 
mean  vertical  advection: 

F.  =  —  +  uQ+f)  +  w— .  (3) 

A  dt  v  JJ  dz 


The  forcing  functions  obtained  for  both  momentum 
and  heat  are  spatially  coherent  (see  section  3b)  and 
qualitatively  consistent  with  similar  diagnostics  presented 
in  other  work  (e.g.,  Bui  et  al.  2009,  their  Fig.  3). 

3)  Other  aspects,  boundary  conditions, 

AND  REGULARIZATION 


To  integrate  the  Sawyer-Eliassen  model,  Eq.  (1)  is 
approximated  with  finite  differences,  as  in  Bui  et  al. 
(2009).  The  radial  and  vertical  grid  spacing  is  2  km  and 
500  m,  respectively.2  A  computational  domain  of  222-km 
radius  and  15-km  height  is  used.  Equation  (1)  is  solved 
numerically  using  the  successive  overrelaxation  scheme 
(SOR)  described  in  Press  et  al.  (1992).  An  overrelaxation 
parameter  of  1.8  is  used  and  the  solution  is  considered 
converged  when  the  absolute  error  of  the  discrete  equiv¬ 
alent  of  Eq.  (1)  is  less  than  the  prescribed  value  of  10-24 
(which  is  typically  10  orders  of  magnitude  smaller  than  the 
magnitude  of  the  forcing). 

The  boundary  conditions  of  Eq.  (1)  at  the  axis  of  rota¬ 
tion  and  at  the  upper  and  lower  boundaries  are  chosen  as 
i[f  =  0  (i.e.,  the  normal  velocity  is  equal  to  zero).  The 
normal  velocity  at  the  outer  radius  (inflow)  is  important  to 
the  conclusions  of  this  study  and  therefore  is  not  prescribed 
(setting  it  to  zero  or  to  the  inflow  value  from  the  RAMS 
integration  would  artificially  constrain  the  solution).  In¬ 
stead,  the  boundary  condition  at  the  outer  radius  is  taken 
as  dif/ldr  =  0  (i.e.,  with  the  vertical  velocity  equal  to  zero). 

Equation  (1)  is  an  elliptic  partial  differential  equation 
given  that  the  discriminant 


D  = 


I  2 


dz 


(xQ 


rd_x 
'  dz 


€x(£+f)  +  c 


(4) 


is  negative  (note  that  the  second  term  on  the  right-hand 
side  is  typically  negative  because  it  is  proportional  to 
dxldz).  At  some  grid  points,  the  balanced  vortices  de¬ 
scribed  in  section  1  (and  the  nonbalanced  vortex  de¬ 
scribed  in  appendix  A)  do  not  satisfy  the  ellipticity 
condition.  The  ellipticity  condition  can  be  violated  when 
d(xC/)dz  is  too  large,  when  [£x(£  +  f)  +  Cd^/dr]  be¬ 
comes  negative,  or  when  dxldz  becomes  positive.  To 
avoid  violation  of  the  ellipticity  condition,  a  regulariza¬ 
tion  procedure  is  implemented,  following  Moller  and 
Shapiro  (2002)  and  Bui  et  al.  (2009). 

As  discussed  in  Moller  and  Shapiro  (2002)  and  in  Bui 
et  al.  (2009),  the  application  of  a  regularization  pro¬ 
cedure  makes  modifications  to  the  stability  parameters 
so  a  convergent  solution  is  obtained  without  changing 
the  basic  vortex  structure.  The  regularization  pro¬ 
cedure  is  as  follows:  In  those  points  where  the  ellipticity 
condition  is  violated,  the  regularization  procedure  re¬ 
duces  the  term  d(xC)/dz  to  0.15  of  its  value.  Then  it 
evaluates  [£x(£  +  f)  +  Cdxldr\  at  every  grid  point.  If 
any  of  the  latter  values  is  negative,  the  magnitude  of  the 
smallest  value  is  multiplied  by  —1.1  and  added  to  all 
grid  points.  Finally,  for  the  balanced  calculations  that 
result  in  negative  vertical  gradients  of  potential  tem¬ 
perature  near  the  top  of  the  boundary  layer  (see  dis¬ 
cussion  below),  at  those  grid  points  where  dxldz  is 
positive  (regions  of  dry,  static  vertical  instability),  the 
procedure  recalculates  the  value  of  potential  temper¬ 
ature  by  assuming  dOldz  =  2Kkm_1.  The  same  ad¬ 
justment  is  applied  to  the  grid  points  above  the  grid 
point  where  the  ellipticity  condition  is  violated,  up  to 
the  height  where  the  adjusted  potential  temperature  is 
less  than  that  of  the  original  vortex.  This  vertical  ex¬ 
tension  of  the  modification  avoids  sharp  gradients  of 
potential  temperature.  As  an  example,  at  80-km  radius, 
this  procedure  modifies  the  temperature  structure 
within  the  layer  between  1.5-  and  4.5-km  heights.  We 
note  that  changing  the  order  of  the  three  criteria  ap¬ 
plied  does  not  change  any  of  the  results  presented.  We 
note  also  that  changing  the  term  d(^C)/dz  to  fractions 
larger  than  0.15  of  its  value  (e.g.,  0.4,  0.6,  or  0.8)  results 
in  local  violation  of  the  ellipticity  criteria,  but  with 
a  solution  that  still  converges  and  renders  a  secondary 
circulation  practically  indistinguishable  from  those 
considered  here. 


3.  Results 


2  A  vertical  grid  spacing  of  300  m  was  also  considered  (not 
shown).  The  higher-resolution  integrations  rendered  solutions  with 
some  small  quantitative  differences  with  respect  to  the  solutions 
shown  below.  However,  none  of  these  differences  altered  in  any 
way  the  conclusions  of  this  manuscript. 


a.  Secondary  eyewall  formation  in  the  RAMS 
integration 

We  begin  this  section  by  reviewing  the  main  evidence 
of  SEF  in  the  RAMS  integration.  Then  we  introduce  the 
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azimuthally  averaged  fields  from  the  RAMS  integra¬ 
tion  and  the  corresponding  calculated  balanced  fields 
used  to  integrate  the  Sawyer-Eliassen  model.  We  de¬ 
scribe  also  details  of  the  Sawyer-Eliassen  integrations 
and  conclude  with  a  discussion  of  departures  from 
balance  dynamics  during  SEF,  as  summarized  in  the 
introduction. 

The  synthetic  storm  in  the  RAMS  integration  evolves 
from  a  weak  mesoscale  convective  vortex  to  a  mature 
hurricane  that  undergoes  a  canonical  eyewall  replacement 
cycle.  Terwey  and  Montgomery  (2008),  Terwey  et  al. 
(2013),  Abarca  and  Montgomery  (2013),  and  Montgomery 
et  al.  (2014)  have  discussed  the  SEF  kinematics  and  basic 
dynamics  of  the  simulation.  The  time  series  of  minimum 
pressure,  radius  of  maximum  winds  (at  about  150-m 
height),  and  maximum  azimuthal-mean  tangential  winds 
(at  1-km  height)  for  the  220-h  numerical  simulation  are 
presented  in  Terwey  and  Montgomery  (2008,  their  Fig.  2). 
The  characteristic  tangential  wind  increase  in  the  SEF 
radial  region  along  with  the  decrease  in  the  tangential 
wind  of  the  primary  eyewall  is  captured  in  the  figure  as 
the  sudden  expansion  of  the  radius  of  maximum  winds 
(from  35-  to  80-km  radius,  just  before  hour  180).  The 
main  characteristics  of  the  azimuthally  averaged  tan¬ 
gential  wind  evolution  and  vertical  velocity  are  pre¬ 
sented  also  in  radius-time  diagrams  in  Terwey  et  al. 
(2013,  their  Fig.  1)  and  in  Abarca  and  Montgomery 
(2013,  their  Fig.  1).  These  diagrams  show  that  the  sim¬ 
ulated  eyewall  replacement  cycle  is  fully  consistent  with 
observations  of  the  phenomena  (e.g.,  Willoughby  et  al. 
1982;  Houze  et  al.  2007;  Bell  et  al.  2012;  Sitkowski  et  al. 
2012)  and  with  other  realistic  numerical  integrations  (e.g., 
Abarca  and  Corbosiero  2011;  Zhou  and  Wang  2011; 
Zhou  et  al.  2011). 

Figure  1  shows  the  radius-height  structure  of  azimuth¬ 
ally  averaged  tangential  wind  velocity  and  secondary  cir¬ 
culation  at  select  times  during  the  RAMS  integration, 
within  200  km  of  the  system-scale  circulation  center.  The 
top  row  corresponds  to  hour  72,  when  the  vortex  was  in¬ 
tensifying  with  a  single  eyewall,  and  the  remaining  rows 
(hours  174, 179, 183,  and  189)  present  the  vortex  with  an 
emerging  and  evolving  secondary  eyewall. 

Figure  1  shows  also  that  the  RAMS  integration  ex¬ 
hibits  the  typical  tangential  wind  structure  and  second¬ 
ary  circulation  of  single-eye  hurricanes  and  a  structure 
fully  consistent  with  the  in  situ  observations  of  hurri¬ 
canes  with  double  eye  walls  (e.g.,  Houze  et  al.  2007,  their 
Fig.  2,  and  Bell  et  al.  2012,  their  Figs.  6  and  8).  Figure  1  is 
similar  to  Abarca  and  Montgomery  (2013),  but  instead 
of  employing  2-h  averages,  we  show  1-h  averaged  fields, 
as  a  simple  way  to  limit  the  impact  of  time  averaging  the 
evolving  flow  during  the  diagnosis  time.  The  right  panels 
allow  for  a  more  direct  comparison  with  observational 


results  (Bell  et  al.  2012,  their  Figs.  6  and  8).  As  empha¬ 
sized  by  Abarca  and  Montgomery  (2013),  the  secondary 
tangential  wind  maximum  emerges  near  the  top,  but 
within  the  boundary  layer.3  We  note  that  the  secondary 
wind  maximum  emerges  in  a  vertically  confined  region, 
with  the  vertical  elongation  of  the  60  m  s_1  isotach  in  Fig. 
lc  spanning  only  about  500  m  in  height.  Such  a  narrow 
vertical  region  with  large  tangential  winds  evolves  to 
occupy  a  broader  vertical  extent,  of  about  2  km  at  hour 
183  and  about  2.5  km  at  hour  189  (Figs.  le,g). 

As  noted  in  Abarca  and  Montgomery  (2013),  Fig.  1 
shows  also  a  much  broader  tangential  wind  field  during 
SEF  (Fig.  lc)  than  during  the  single -eyewall  phase  of 
storm  intensification  (Fig.  la).  Abarca  and  Montgomery 
(2013)  showed  that  this  structure  change  is  a  progressive 
one  (see  their  Fig.  la),  consistent  with  the  ideas  presented 
by  Huang  et  al.  (2012,  their  Figs,  lc-f)  and  Wu  et  al. 
(2012,  their  Fig.  6a).  Finally,  the  figure  shows  how  the 
forming  secondary  eyewall  is  characterized  by  a  broad 
radial  region  of  ascending  motion  (Fig.  Id)  that  becomes 
radially  narrower  with  time  (Figs.  lf,h)  as  the  secondary 
eyewall  forms.  This  progressive  radial  focusing  of  the 
azimuthally  averaged  upward  motion  is  also  visible  in 
Fig.  lb  of  Abarca  and  Montgomery  (2013)  and  seems  to 
be  a  characteristic  signature  of  SEF. 

b.  Balanced  temperature  and  wind  fields 

We  now  present  and  discuss  balanced  fields  calcu¬ 
lated  from  azimuthally  averaged  and  time-averaged 
RAMS  fields,  as  captured  during  SEF.  As  a  represen¬ 
tative  time,  we  show  the  results  centered  at  hour  174. 
However,  the  described  processes  are  consistent  with 
findings  from  other  analysis  times  (e.g.,  hours  162, 166, 
170;  not  shown).  Figure  2  shows  the  azimuthally  av¬ 
eraged  tangential  wind  field  at  174  h  (as  in  Fig.  lc)  and 
its  corresponding  balanced  potential  temperature 
field.  The  latter  is  calculated  following  (Smith  2006),  as 
detailed  in  section  2c(l),  using  the  tangential  wind  field 
in  Fig.  2a.  As  expected  from  the  thermal  wind  equa¬ 
tion,  the  balanced  potential  temperature  field  exhibits 
a  warm-core  structure  in  most  of  the  troposphere, 
above  the  height  of  maximum  winds.  It  shows  also 
a  cold-core  structure  at  very  low  levels,  in  the  bound¬ 
ary  layer,  in  the  region  where  tangential  winds  increase 
with  height. 

In  the  balanced  field,  the  near-surface  potential 
temperature  at  around  30-km  radius  (below  the  loca¬ 
tion  of  the  wind  maximum)  reaches  a  value  lower  than 
-30°C.  The  unrealistic  cold  anomaly  was  observed  also 


3  See  the  introduction  for  our  working  definition  for  the 
boundary  layer. 
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Fig.  1.  Radius-height  sections  of  azimuthally  averaged  (left)  tangential  (contours)  and  (right)  vertical  and  radial  (vectors)  components 
of  the  storm-relative  wind.  Five  different  integration  times  are  shown,  corresponding  to  hours  (a),(b)  72,  (c),(d)  174,  (e),(f)  179,  (g),(h)  183, 
and  (i),(j)  189.  Hours  174  and  179  are  referred  to  as  hours  18  and  23  in  Terwey  and  Montgomery  (2008),  Terwey  et  al.  (2013),  and  Abarca 
and  Montgomery  (2013).  Note  that  (a),(b)  correspond  to  about  100  simulation  hours  before  the  rest  of  the  panels.  Such  an  early  time  is 
presented  to  contrast  the  wind  field  during  a  time  of  single-eyewall  intensification  and  during  SEF  (see  text).  Note  also  that  (e)  corre¬ 
sponds  to  Fig.  2f  in  Abarca  and  Montgomery  (2013),  except  here  a  snapshot  (instead  of  2-h  average)  is  shown.  Contours  are  shown  every 
5  m  s_1  and  the  reference  vector  magnitude  is  indicated  below  (j).  The  vectors  are  curved  so  they  are  tangential  to  the  represented  flow  at 
each  point.  Note  that  (c)-(f)  are  also  published  as  in  the  top  panels  of  Fig.  2  of  Montgomery  et  al.  (2014). 


in  Bui  et  al.  (2009,  their  Fig.  2d),  with  their  potential 
temperature  field  reaching  values  lower  than  7°C  at  the 
lowest  part  of  their  domain.  For  the  case  studied 
herein,  the  cold  anomaly  is  exacerbated  both  by  the 
intensity  of  the  simulated  storm  [maximum  azimuthally 
averaged  tangential  winds  surpassing  70 ms-1,  com¬ 
pared  with  about  50 ms-1  in  Bui  et  al.  (2009)]  and 


by  the  broad  radial  structure  of  the  tangential  wind 
field.4  The  unrealistic  cold  anomaly  is  found  also  by 


4  Observe  that  in  Bui  et  al.  (2009,  their  Fig.  2b)  the  50m  s_1 
isotach  spans  from  about  35-  to  about  45-km  radius,  while  the  same 
isotach  roughly  spans  from  22-  to  120-km  radius  in  Fig.  2a. 
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Corresponding  balanced  potential  temperature  [K] 


Fig.  2.  Radius-height  sections  of  (a)  azimuthally  averaged  tan¬ 
gential  component  of  the  storm-relative  wind  (contours  every 
5  m  s_1)  for  hour  174  (as  in  Fig.  lb)  and  (b)  corresponding  balanced 
temperature  field  (contours  every  10  K)  calculated  as  in  Smith 
(2006;  see  text  for  details). 

Rozoff  et  al.  (2012).  Given  such  an  unrealistic  cold 
anomaly,  Rozoff  et  al.  (2012)  chose  to  hold  the  tan¬ 
gential  winds  constant  below  2  km. 

The  unrealistic  cold  anomaly  in  the  near-surface  po¬ 
tential  temperature  field  in  balance  with  the  tangential 
wind  is  in  stark  contrast  with  the  RAMS  azimuthally 
averaged  potential  temperature  in  that  region  (32°C; 
Fig.  3a).  Besides  the  large  cold  anomaly  at  the  surface, 
the  balanced  temperature  field  derived  using  Smith 
(2006)  exhibits  static  instability  in  the  inner  core  as 
a  result  of  the  change  of  sign  of  the  vertical  wind  shear 
near  the  top  of  the  boundary  layer  and  its  implications 
through  the  thermal  wind  equation.  In  contrast,  strong 
changes  in  the  vertical  stability  are  not  a  characteristic  of 
the  azimuthally  averaged  temperature  field  (Fig.  3a).  It 
is  a  useful  reminder  to  note  that  the  unrealistic  cold 
anomaly  and  vertical  stability  profiles  result  from  the 
assumption  that  the  potential  temperature  field  is  in 
balance  with  the  prescribed  tangential  wind  in  the 
boundary  layer.  Since  we  aim  to  assess  departures  from 
balance,  we  will  adhere  fully  to  this  assumption.  While 
the  magnitude  of  the  temperature  difference  can  be 
thought  as  a  measure  of  the  inappropriateness  of  the  as¬ 
sumption  of  balance  in  the  boundary  layer,  it  permits  the 
representation  of  a  truly  balanced  vortex  in  the  Sawyer- 
Eliassen  inversion. 


Azimuthal  average  potential  temperature  [K] 


Fig.  3.  Radius-height  sections  of  (a)  the  azimuthally  averaged 
potential  temperature  for  hour  174  and  (b)  the  corresponding  bal¬ 
anced  tangential  wind  velocity,  calculated  as  in  Holton  (2004;  see 
text).  Contours  are  shown  every  5ms_1  and  5K,  except  in  (b)  for 
radii  larger  than  80  km  and  heights  lower  than  9  km,  where  contours 
are  shown  every  0.5 ms-1.  Note  that  (b)  is  also  published  as  in  the 
left-center  panel  of  Fig.  2  of  Montgomery  et  al.  (2014). 

Figure  3  shows  the  azimuthally  averaged  potential 
temperature  field  from  the  RAMS  integration,  along 
with  the  tangential  wind  field  in  balance  with  the  RAMS 
azimuthally  averaged  density  and  pressure  fields  (not 
shown).  Here,  the  tangential  wind  field  is  computed 
from  the  radial  pressure  gradient  force,  neglecting  the 
radial-flow  component,  as  is  customary  in  the  traditional 
balance  formulation,  as  described  in  section  2c(l).  In  the 
lower  troposphere,  there  are  acute  differences  between 
the  azimuthally  averaged  tangential  wind  and  the  anal¬ 
ogous  balanced  field.  The  tangential  velocity  maximum 
changes  from  more  than  80  m  s-1  in  the  average  to  about 
70 ms-1  in  the  balanced  field  (the  balanced  tangential 
wind  maximum  is  about  80%  of  the  magnitude  of  the 
azimuthally  averaged  wind  maximum),  and  it  is  located 
about  6  km  radially  outward  and  1.5  km  above  that  of  the 
azimuthally  averaged  maximum.  Despite  these  differences 
in  the  velocity  maximum,  the  largest  discrepancies  be¬ 
tween  the  azimuthal  average  and  the  balanced  tangen¬ 
tial  wind  field  occur  radially  outward  of  the  radius  of 
maximum  winds,  in  the  SEF  region.  The  balanced  tan¬ 
gential  velocity  field  Fig.  3b  exhibits  three  discernible 
tangential  wind  maxima  outside  of  the  primary  eyewall. 
This  is  in  contrast  with  the  RAMS  data  (Fig.  2a),  which 
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have  no  tangential  wind  maximum  beyond  the  primary 
eyewall  at  hour  174  and  possess  just  one  secondary 
tangential  wind  maximum  at  about  5h  later  (Fig.  lg). 
The  three  maxima  in  the  balanced  tangential  wind  field 
are  located  in  the  lowest  2.5  km  of  the  model  atmosphere, 
where  the  differences  between  the  azimuthal  average  and 
the  corresponding  balanced  fields  are  not  small. 

Just  like  the  stark  contrast  between  the  RAMS  po¬ 
tential  temperature  (Fig.  3a)  and  the  balanced  potential 
temperature  field  (Fig.  2b),  the  substantial  differences 
between  the  RAMS  tangential  wind  field  (Fig.  2a)  and 
the  computed  balanced  tangential  wind  field  (Fig.  2b)  are 
quantitative  indications  of  the  inappropriateness  of  the 
assumption  that  hurricanes  are  in  balance  at  all  heights. 
The  differences  are  largest  in  the  boundary  layer  where 
the  tangential  wind  maximum  of  the  secondary  eyewall 
emerges. 

c.  Forcing  fields 

To  integrate  the  Sawyer-Eliassen  model,  the  diabatic 
heating  rate  and  tangential  momentum  forcing  functions 
are  required.  These  functions  are  computed  from  azi¬ 
muthal  and  1-h  averages  of  the  RAMS  integration  as 
described  in  section  2c(2).  Figure  4  shows  the  forcing 
functions  used  in  this  study.  Consistent  with  the  ob¬ 
served  secondary  circulation  (Fig.  Id),  Fig.  4a  shows  that 
the  azimuthally  averaged  diabatic  heating  rate  exhibits 
a  well-defined  maximum  confined  largely  to  a  vertical 
column  within  an  annular  region  about  15  km  wide, 
centered  around  40  km  and  sloping  outward.  The  mean 
heating  rate  maximum  exhibits  a  magnitude  of  about 
70  K  h_1.  This  value  is  substantially  higher  than  those  pre¬ 
viously  reported  [e.g.,  about  20Kh_1  in  Bui  et  al.  (2009)]. 
Both  the  strong  vortex  (with  tangential  wind  maximum 
larger  than  80 ms-1)  and  a  highly  symmetric  structure 
[Fig.  4  of  Terwey  and  Montgomery  (2008)]  are  factors 
contributing  to  this  large  value.  As  the  primary  eyewall 
begins  to  weaken  around  this  time,  the  azimuthally  aver¬ 
aged  diabatic  heating  rate  decreases  substantially  reach¬ 
ing  values  of  60Kh_1  at  hour  179.  Figure  4a  shows  also 
that  as  the  vortex  undergoes  SEF,  the  heat  forcing  exhibits 
rates  that  exceed  10  K  h_1  between  70-  and  150-km  radius. 
Figure  4b  shows  a  sink  of  tangential  momentum  in 
a  shallow  layer  above  the  surface,  attributable  to  the  ef¬ 
fects  of  surface  friction  that,  by  itself,  is  acting  persistently 
to  decelerate  the  tangential  wind.  Not  surprisingly,  the 
maximum  deceleration  is  found  in  the  region  with  the 
largest  tangential  wind  speeds.  The  deceleration  reaches 
50m s1  h1  in  the  region  of  the  primary  eyewall  and  sur¬ 
passes  10  m  s_1  h_1  at  around  100-km  radius.  In  the  middle 
and  upper  troposphere,  Fig.  4b  shows  that  there  is  a  mo¬ 
mentum  source  radially  inwards  of  the  primary  eyewall 
updraft  (Fig.  Id)  and  a  momentum  sink  radially  outwards 


Radius  [km] 


Fig.  4.  Radius-height  sections  of  azimuthally  averaged  forcing 
functions  diagnosed  at  hour  174  from  the  three-dimensional 
RAMS  calculation:  (a)  the  heat  source  (contours  every  10  Kh1) 
and  (b)  the  momentum  source  (contours  every  10ms_1h_1).  Neg¬ 
ative  contours  are  dashed. 


of  the  same  updraft.  These  sources  and  sinks  are  related 
to  eddy  momentum  fluxes  in  the  eyewall  region.  While 
the  magnitude  of  this  forcing  is  large  locally  (about 
20ms_1h_1),  it  is  relatively  small  compared  to  the 
forcing  at  lower  levels,  near  the  surface.  As  shown 
below,  the  effect  of  this  eddy  momentum  forcing  on  the 
overall  solution  for  the  mean  transverse  circulation  is 
relatively  small. 

d.  RAMS  versus  Sawyer-Eliassen  tangential  wind 
tendency 

Figure  5  compares  the  azimuthally  averaged  tan¬ 
gential  wind  tendency  and  the  radial,  vertical,  and 
tangential  wind  velocity  of  the  RAMS  integration  at 
hour  174  (left  column),  with  those  resulting  from  the 
Sawyer-Eliassen  model,  integrated  as  described  in 
section  3b,  following  Smith  (2006)  (center  column)  and 
Holton  (2004)  (right  column).  The  tangential  wind 
tendency  (lowest  two  rows)  is  computed  as  indicated  in 
the  figure,  as 

^=-i7(£+/)-W^  (5) 

ot  dZ 

for  the  third  row,  and  as 
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Fig.  5.  Vertical  cross  sections  of  the  (a)-(c)  radial  and  (d)-(f)  vertical  wind  components  and  (g)-(l)  tangential  wind  tendencies  (defined 
in  the  equations  on  the  left  side  of  the  figure),  for  (left)  the  azimuthally  averaged  fields  of  the  RAMS  simulation  at  hour  174  and  for  the 
Sawyer-Eliassen  model  integrated  with  balanced  fields  as  in  (center)  Smith  (2006)  and  (right)  Holton  (2004).  See  text  for  further  details 
and  definition  of  symbols  in  the  equations.  Contour  intervals  are  indicated  above  each  panel. 
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for  the  bottom  row.  In  these  two  equations,  the  symbols 
are  as  defined  before.  In  the  case  of  the  Sawyer-Eliassen 
model  (center  and  right  columns  in  Fig.  5),  we  compute 
the  tangential  wind  tendencies  in  Eqs.  (5)  and  (6)  using 
the  balanced  secondary  circulation  deduced  from  Eq. 
(1)  instead  of  the  azimuthal  averages  u  and  w  from  the 
RAMS  solution.  Note  that  in  the  center  and  right  col¬ 
umns  in  Fig.  5,  the  difference  between  the  third  and 
bottom  row  lies  with  the  inclusion  of  the  azimuthally 
averaged  tangential  momentum  sink. 

Figure  5  shows  that  both  the  RAMS  and  the  Sawyer- 
Eliassen  results  exhibit  a  secondary  circulation  typical  of 
a  hurricane  undergoing  SEF,  with  inflow  in  the  bound¬ 
ary  layer  and  in  the  midtroposphere,  and  outflow  in  the 
upper  troposphere  and  immediately  above  the  strong 
inflow  in  the  boundary  layer.  Figure  5  shows  also  that 
while  the  Sawyer-Eliassen  inversions  capture  the  over¬ 
all  secondary  circulation  structure  of  the  RAMS  simu¬ 
lation,  there  are  substantial  differences  between  the 
RAMS  and  Sawyer-Eliassen  results.  As  an  example,  the 
largest  magnitude  of  vertical  motion  in  the  three  com¬ 
putations  is  observed  within  the  primary  eyewall  (around 
40-km  radius  in  the  three  cases;  see  Figs.  5d-f),  with  the 
balance  computation  following  Smith  (2006)  exhibiting 
the  largest  values  (surpassing  4.5  m  s_1)  and  the  smallest 
(<2.5ms_1)  in  the  balance  computation  following 


Holton  (2004).  The  larger  upward  vertical  motion  in  the 
balance  calculation  following  Smith  (2006)  exhibits  also 
larger  adjacent  subsidence  (>1.0ms_1).  Such  strong 
subsidence  is  not  observed  in  the  azimuthally  averaged 
RAMS  field  or  in  the  balance  calculation  following 
Holton  (2004). 

In  the  region  of  SEF  (between  75-  and  150-km  radius) 
the  RAMS  maximum  azimuthally  averaged  boundary 
layer  inflow  surpasses  values  of  16 ms-1  (Fig.  5a).  In  con¬ 
trast,  the  Sawyer-Eliassen  inversions  exhibit  much  smaller 
mean-inflow  velocities  that  barely  surpass  4ms”1  (Fig.  5b) 
for  the  case  of  the  Smith  (2006)  balance  methodology 
and  3 ms-1  (Fig.  5c)  for  the  case  of  the  Holton  (2004) 
balance  methodology.  In  the  same  radial  region,  the  ver¬ 
tical  motion  in  the  RAMS  data  exhibits  an  azimuthal 
average  with  magnitude  larger  than  0.5 ms-1  (Fig.  5d). 
The  balance  calculation  computed  using  the  Smith 
(2006)  pathway  exhibits  also  some  localized  upward  verti¬ 
cal  motion  in  the  region  of  SEF  that  surpasses  the  0.5  ms-1 
level  (see  contours  at  about  86-km  radius  and  3-km 
height  in  Fig.  5e),  unlike  the  balance  calculation  com¬ 
puted  using  the  Holton  (2004)  pathway  that  does  not 
reach  such  values. 

Figure  5  shows  that  the  differences  between  the  azi¬ 
muthally  averaged  radial  and  vertical  RAMS  velocity 
fields  and  those  of  the  Sawyer-Eliassen  model  in¬ 
tegrations  are  not  of  overall  structure  but  of  magnitude. 
The  inflow  in  the  boundary  layer  is  substantially  under¬ 
estimated  in  the  Sawyer-Eliassen  calculations,  by  up  to 
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about  80%  in  the  radial  region  of  SEF.  The  figure  shows 
also  that  quantitative  differences  in  the  secondary  circu¬ 
lation  translate  to  quantitative  and  qualitative  differences 
in  tangential  wind  tendency,  with  none  of  the  Sawyer- 
Eliassen  results  being  able  to  capture  even  the  right  sign 
of  the  tangential  wind  tendency  in  the  region  of  SEF  (see 
appendix  B  for  yet  another  pathway  to  characterize  the 
balanced  fields  consistent  with  these  results).  In  the  radial 
region  of  SEF  (between  about  75-  and  150-km  radius), 
the  tangential  wind  tendency,  computed  from  the  sec¬ 
ondary  circulation  using  the  Sawyer-Eliassen  integration 
renders  only  a  modest  tangential  wind  tendency  (com¬ 
posed  of  a  weak  radial  absolute  vorticity  flux  and  vertical 
advection  of  tangential  momentum;  see  Figs.  5h,i).  When 
the  azimuthally  averaged  momentum  sink  [Eq.  (3)]  is 
added  to  the  foregoing  tangential  wind  tendency  com¬ 
putation,  based  on  the  secondary  circulation  of  the  bal¬ 
ance  integration,  no  SEF  spinup  tendency  is  predicted. 
The  quantitative  differences  in  inflow  result  in  qualitative 
differences  in  tangential  wind  tendency,  since  ultimately 
the  tangential  wind  spinup  (in  the  azimuthally  averaged 
sense)  is  due  largely  to  the  radial  flux  of  absolute  vorticity 
[see  Abarca  and  Montgomery  (2013,  section  5c,  their  tan¬ 
gential  velocity  tendency  diagnosis)].  Figure  5  shows  that 
the  underestimation  of  the  inflow  magnitude,  and  conse¬ 
quently  the  underestimation  of  the  absolute  vorticity  flux, 
when  combined  with  the  vertical  advection  of  tangential 
momentum,  results  in  spindown  in  the  radial  region  where 
SEF  is  taking  place. 

The  results  presented  here  suggest  that  during  the 
simulated  SEF  event,  the  magnitude  of  the  radial  inflow 
needed  to  maintain  the  specified  vortices  in  axisymmetric 
balance  (as  defined  herein)  is  much  smaller  than  the  in¬ 
flow  generated  by  latent  heating  and  surface  friction  in  the 
full-physics  simulation.  The  full-physics  simulation  in¬ 
cludes  strong  horizontal  advective  dynamics  that  arise  in 
the  region  spanning  the  transition  layer  between  the  sur¬ 
face  layer  and  the  nearly  inviscid  bulk  vortex  circulation. 
Because  of  the  strongly  coupled  nature  of  the  horizontal 
momentum  equations  in  the  composite  boundary  layer 
region,  the  frictional  force  in  the  boundary  layer  results  in 
subgradient  winds  near  the  surface,  where  the  radial 
pressure  gradient  force  is  larger  than  the  sum  of  the 
centrifugal  and  Coriolis  forces.  The  subgradient  mo¬ 
mentum  balance  near  the  surface  results  in  an  increase 
of  the  radial  inflow  and  a  correspondingly  large  radial 
vorticity  flux,  the  net  result  of  which  can  be  sufficient  to 
generate  supergradient  flow  and  tangential  wind 
spinup  farther  up  within  the  boundary  layer.  These 
results  suggest  that  balance  dynamics  (e.g.,  Nong  and 
Emanuel  2003;  Rozoff  et  al.  2012)  are  generally  an  in¬ 
sufficient  theoretical  framework  to  explain  secondary 
eyewall  formation. 


44 


I 


~T~ 

50 


Radius  [km] 


Fig.  Al.  Radius-height  sections  of  the  discriminant  D.  Gray 
shading  indicates  regions  where  the  ellipticity  condition  is  violated 
(D  <  0)  for  the  integrations  performed  by  calculating  the  balance 
as  in  (a)  Smith  (2006)  and  (b)  Holton  (2004). 


4.  Conclusions 

In  this  study  we  examine  departures  from  axisym¬ 
metric  balance  dynamics  during  secondary  eyewall  for¬ 
mation,  in  a  case  of  a  three-dimensional  full-physics 
integration  of  a  tropical  cyclone.  The  mesoscale  simu¬ 
lation  is  convection  permitting,  it  is  integrated  in  an 
idealized  quiescent  environment  from  a  weak  mesoscale 
vortex,  and  it  evolves  to  develop  a  canonical  eyewall 
replacement  cycle.  This  study  focuses  on  the  time  of  the 
simulation  during  the  formation  of  the  secondary  eye- 
wall,  a  few  hours  before  the  secondary  tangential  wind 
maxima  are  established. 

Departures  from  axisymmetric  balance  dynamics  are 
quantified  by  comparing  the  secondary  circulation  and 
corresponding  tangential  wind  tendency  of  the  mesoscale 
simulation  with  those  diagnosed  after  constraining  the 
azimuthally  averaged  flow  to  balance  dynamics.  We  refer 
to  balance  dynamics  as  those  processes  that  act  to  main¬ 
tain  axisymmetric  thermal  wind  balance.  In  contrast,  we 
refer  to  unbalanced  dynamics  as  those  processes  that  are 
related  to  the  axisymmetric  portion  of  the  hurricane  flow 
that  are  not  constrained  by  the  thermal  wind  equation.  In 
our  view,  unbalanced  dynamics  include  the  strong  hori¬ 
zontal  advective  dynamics  that  arise  in  hurricanes  on 
account  of  surface  friction  and  therefore  are  confined 
primarily  to  the  low-level  region  spanning  the  transition 
layer  between  the  surface  layer  and  the  nearly  inviscid 


3734 


JOURNAL  OF  THE  ATMOSPHERIC  SCIENCES 


Volume  71 


dv 

di 


XI 

SjO  e 

’5 
®  3J 


dv  r  ■ 


Radius  [km] 


Fig.  Bl.  Radius-height  sections  of  (a)  the  discriminant,  (b)  radial  and  (c)  vertical  wind 
components,  and  (d)  and  (e)  tangential  wind  tendencies  (defined  in  the  equations  on  the  left 
side  of  the  figure),  for  the  Sawyer-Eliassen  model  integrated  with  azimuthally  averaged  fields 
(not  in  balance)  from  the  RAMS  simulation  at  hour  174.  See  text  for  details.  Contour  intervals 
are  indicated  above  each  panel. 


bulk  vortex  circulation.  We  use  the  Sawyer-Eliassen 
balance  equation,  commonly  used  as  a  mean-flow  model 
in  the  limit  of  weak  departures  from  circular  flow,  to 
diagnose  the  secondary  circulation  constrained  by  axi- 
symmetric  balance  dynamics.  To  integrate  the  Sawyer- 
Eliassen  equation  we  characterize  the  background  vortex 
and  diagnose  the  heat  and  momentum  forcing  from  the 
mesoscale  numerical  integration. 

Three  different  ways  of  characterizing  the  basic-state 
vortex  in  the  Sawyer-Eliassen  integration  are  studied.  In 


one  method,  we  use  the  azimuthally  averaged  tangential 
wind  field  and  compute  its  corresponding  balanced  ther¬ 
modynamical  fields,  using  the  thermal  wind  equation.  In 
a  second  method,  we  use  the  azimuthally  averaged  pres¬ 
sure  and  density  fields  from  the  mesoscale  integration  and 
infer  the  corresponding  balanced  tangential  velocity  field. 
Finally,  as  a  third  method,  we  characterize  the  vortex 
using  only  azimuthally  averaged  fields  from  the  mesoscale 
integration.  Azimuthal  averages  and  their  analogous  es¬ 
timated  balanced  fields  exhibit  many  commonalities 
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throughout  the  troposphere.  However,  these  fields  exhibit 
large  differences  in  the  boundary  layer.  Such  differences 
are  an  alternative  indication  of  the  inappropriateness  of 
the  assumption  of  axisymmetric  balance  dynamics  (or 
weak  departures  therefrom)  in  the  boundary  layer. 

The  most  significant  finding  of  this  work  is  that  axi¬ 
symmetric  balance  dynamics  does  not  capture  the  tan¬ 
gential  wind  spinup  during  secondary  eyewall  formation. 
While  balance  dynamics  does  capture  the  overall  structure 
of  the  secondary  circulation  of  the  full-physics  simula¬ 
tion,  it  exhibits  quantitative  differences  that  include 
a  large  underestimation  of  inflow  in  the  boundary  layer, 
where  the  secondary  wind  maximum  emerges.  The 
quantitative  underestimation  of  inflow  results  in  striking 
qualitative  differences  in  the  tangential  wind  tendency, 
and  balance  dynamics  are  found  to  be  unable  to  capture, 
by  themselves,  the  secondary  eyewall  spinup.  These 
findings  are  robust  to  the  pathway  of  balanced  vortex 
characterization  for  the  three  different  methodologies 
used  in  this  study. 

The  results  herein  offer  support  to  previous  studies 
that  point  to  the  importance  of  unbalanced  boundary 
layer  dynamics  in  secondary  eyewall  formation.  These 
results  suggest  that  balance  dynamics  are  generally  an 
insufficient  theoretical  framework  to  explain  secondary 
eyewall  formation.  Ongoing  analysis  that  parallels  the 
work  presented  here  is  being  conducted  on  other  high- 
resolution  numerical  simulations  carried  out  using  the 
Weather  Research  and  Forecasting  (WRF)  Model.  The 
results  of  this  work  will  be  presented  in  due  course. 
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APPENDIX  A 

Violation  of  the  Ellipticity  Condition 

The  balanced  vortex  described  in  Figs.  2  and  3,  forced 
with  the  functions  in  Fig.  4  represents  an  elliptic  problem 
in  most  points  in  the  domain.  Figure  A1  shows  the  areas 
where  the  ellipticity  condition  is  met  (white  area)  and 
where  it  is  violated  (gray  area),  for  the  balanced  vortex 


defined  as  in  Smith  (2006)  (Fig.  Ala)  and  as  in  Holton 
(2004)  (Fig.  Alb).  The  two  approaches  to  define  balance 
result  in  local  violations  of  the  ellipticity  condition  in  the 
upper  part  of  the  domain.  These  violations  are  related  to 
[£x(£  +  f)  +  Cdx/dr]  becoming  negative  and  mostly 
occur  in  the  region  of  outflow.  The  two  approaches  ex¬ 
hibit  also  some  violations  of  the  ellipticity  condition  at 
large  radius,  near  the  lower  boundary,  related  to  a  large 
d(xC)/dz.  Besides  these  regions  of  violation  of  the  el¬ 
lipticity  condition,  the  computation  based  on  Smith 
(2006)  results  also  in  violations  in  the  lower  troposphere 
near  the  storm  center  (between  1-  and  3-km  heights  and 
for  radii  less  than  80  km)  related  to  d^/dz  becoming 
positive.  The  latter  violations  are  related  to  the  vertical 
instability  of  the  computed  balanced  temperature  field. 
With  the  regularization  procedure  described  in  section 
3,  the  ellipticity  condition  is  met  throughout  the  domain 
and  the  solution  of  Eq.  (1)  converges  (absolute  error 
smaller  than  10-24)  after  a  few  thousand  iterations. 

APPENDIX  B 

Alternative  Vortex  Representation 

In  this  appendix  we  compare  the  secondary  circula¬ 
tion  and  the  corresponding  tangential  wind  tendencies 
of  the  RAMS  and  Sawyer-Eliassen  models,  using  al¬ 
ternative  vortex  representations.  To  define  the  basic- 
state  vortex,  the  Sawyer-Eliassen  equation  is  integrated 
using  azimuthally  averaged  dynamical  and  thermody¬ 
namical  fields  (Figs.  2a  and  3a). B1  The  forcing  functions 
employed  are  those  in  Fig.  4.  With  this  configuration,  the 
diagnostic  balance  equation  represents  an  elliptic 
problem  in  most  points  in  the  domain.  Figure  Bla  shows 
the  areas  where  the  ellipticity  condition  is  met  (white 
area)  and  where  it  is  violated  (gray  area).  Figure  Bla 


B1  Another  way  of  representing  the  vortex,  essentially  combining 
the  two  methods  in  the  main  body  of  this  work,  has  been  considered 
following  a  suggestion  by  D.  Nolan  (2014,  personal  communica¬ 
tion).  This  procedure  consists  of  taking  the  vortex  from  the  second 
method  (Holton  2004;  see  main  text)  and  computing  the  thermo¬ 
dynamical  fields  that  would  be  in  thermal  wind  balance  with  its 
wind  field  (through  the  first  method;  Smith  2006;  see  main  text). 
Strictly  speaking,  this  second  iteration  is  required  to  ensure  con¬ 
sistency  between  the  mass  and  wind  fields  using  the  Holton  (2004) 
methodology.  With  the  second  iteration,  the  corrections  to  the 
temperature  field  are  relatively  small  through  most  of  the  vortex, 
with  the  largest  corrections  located  in  the  upper  troposphere  within 
approximately  50-km  radius  from  the  center  of  the  storm  (not 
shown).  With  the  second  iteration,  the  results  from  the  Sawyer- 
Eliassen  inversion  and  the  derived  tangential  wind  tendency  (not 
shown)  are  qualitatively  and  quantitatively  very  similar  to  that  of 
the  second  method  (Holton  2004;  see  main  text). 


3736 


JOURNAL  OF  THE  ATMOSPHERIC  SCIENCES 


Volume  71 


shows  that  the  ellipticity  condition  is  violated  locally,  at 
the  lower  part  of  the  domain.  This  is  related  to  a  large 
d(xC)ldz,  as  described  in  section  3.  With  the  regulari¬ 
zation  procedure  as  described  in  section  2c(3),  the  el¬ 
lipticity  condition  is  met  throughout  the  domain  and  the 
solution  of  Eq.  (1)  converges  (error  smaller  than  10-24) 
after  a  few  thousand  iterations. 

Figures  Bib  and  Blc  show  the  secondary  circulation 
resulting  from  the  Sawyer-Eliassen  integration.  These 
figures  show  that  the  spatial  structure  of  the  secondary 
circulation  of  the  corresponding  full-physics  computa¬ 
tion  is  again  broadly  captured.  The  balanced  secondary 
circulation  even  exhibits  a  near-surface,  secondary  in¬ 
flow  maxima  that  surpasses  4 ms-1  around  90-km  ra¬ 
dius.  However,  the  inflow  values  in  the  SEF  region  are 
much  smaller  than  those  in  the  full-physics  integration 
(cf.  Fig.  5a). 

Figure  Bid  shows  the  tangential  wind  tendency 
computed  from  the  Sawyer-Eliassen  solution,  using 
only  the  radial  absolute  vorticity  flux  and  the  vertical 
advection  of  tangential  momentum.  The  figure  shows 
modest  positive  tangential  wind  tendency  in  the  radial 
region  of  SEF  (between  about  75-  and  150-km  radius). 
However,  Fig.  Ble  shows  that  when  the  azimuthally 
averaged  momentum  sink  [Eq.  (3)]  is  added  to  the 
tangential  wind  tendency,  no  spinup  tendency  is  ob¬ 
tained  in  the  radial  region  of  SEF.  That  is,  quantitative 
differences  in  the  secondary  circulation  between  the 
full-physics  model  (Fig.  5,  left  column)  and  the  balance 
diagnostic  equation  represent  striking  qualitative  dif¬ 
ferences  in  storm  spinup,  as  in  the  results  presented  in 
section  3.  The  foregoing  results  support  the  view  that 
axisymmetric  balance  dynamics  as  defined  herein  does 
not  capture  the  tangential  wind  spinup  during  second¬ 
ary  eyewall  formation. 
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